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Abstract 

We recently presented a new method for the evaluation of one-loop amplitude of arbitrary 
scattering processes, in which the reduction to scalar integrals is performed at the integrand 
level. In this talk, we review the main features of the method and briefly summarize the 
results of the first calculations performed using it. 
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1 Introduction and General Motivations 



The experimental programs of LHC require high precision predictions for multi-particle pro- 
cesses. At the tree level, the introduction of efficient recursive algorithms pQ improved the 
theoretical description of such processes. However, the current need for precision goes beyond 
tree order. The search and the interpretation of new physics requires a precise understand- 
ing of the Standard Model. We need accurate predictions and reliable error estimates. For 
future experiments, starting with LHC, all analyses will require at least next-to- leading order 
calculations (NLO). 

In the last few years, several groups have been working on the problem of constructing 
efficient and automatized methods for the computation of one-loop corrections for multi-particle 
processes. 

The standard reference for any one-loop calculation is the fundamental work of 't Hooft 
and Veltman[2], and Passarino and Veltman[3], which already contains many of the ingredients 
needed to accomplish such calculations. However, after almost three decades, only few one- loop 
calculations involving more than five particles have been completed [4]. 

The difficulties arising in this kind of calculations are well known: on the one hand the 
presence of a very large number of Feynman diagrams, on the other the appearance of numerical 
instabilities (i.e Gram determinants), that should be cured or avoided. 

Many different interesting techniques have been proposed for NLO calculations: numerical 
methods [5], in which tensor integrals are directly computed numerically; semi-numerical meth- 
ods [6], in which a reduction to a basis of known integrals is performed, dealing carefully with 
spurious singularities; analytic approaches [7], that make use of unitarity cuts to build NLO 
amplitudes by gluing on-shell tree amplitudes. Some of these techniques require additional ra- 
tional terms to be computed separately [8]. Recent complete reviews of existing methods can 
be found, for example, in Refs. [91 [10] . 

The main purpose of this talk is to illustrate the new method of reduction at the integrand 
level (OPP reduction) that we developed during the last year [11|,|12]. The method benefits from 
previous work of Pittau and del Aguila [13]. In Section 2, we introduce the main features of the 
method. In Section 3, we deal with the computation of rational terms. Finally, in Section 4, 
we summarize the applications of the method that have already been implemented. 



2 OPP Reduction 

Any m-point one-loop amplitude can be written, before integration, as 
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with Di = (q + Pi) 2 — mf. The bar denotes objects living in n = 4 + e dimensions and a 
tilde objects of dimension e. Physical external momenta pi are 4-dimensional objects, while 
the integration momentum q is in general n-dimensional. Following this notation, we have 
f = q 2 + q 2 and A = A + q 2 ■ 



Assuming for the moment that the numerator N(q) is fully 4-dimensional, we can rewrite 
it at the integrand level in terms of D{ as 
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The quantities <i(&o*l*223) are the coefficients of 4-point scalar functions with denominators 
labeled by io, i±, 12, and 13. In the same way, we call c{ioi\i2), K^oh), «(^o) the coefficients of 
all possible 3-point, 2-point and 1-point scalar functions, respectively. 

The quantities d, c, b, a are what we define "spurious" terms, i.e. terms that are present in 
the decomposition at the integrand level, but they will vanish upon integration. These terms 
still depend on the integration momentum q. 

In order to perform the reduction, we should first find a general explicit expression for the 
spurious term. Following Ref. [13], we rewrite all q's in N(q) by means of a basis of massless 
vectors if 

4 
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where 

fcl = l\ + ail 2 , k 2 = £2 + a 2 £i , h = pi - p 

and 
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The coefficients G{ either reconstruct denominators Di or vanish upon integration. In the first 
case, they give rise to d, c, b, a coefficients, while in the second one they form the spurious d, 
c, b, a coefficients. 

As an example, we explicitly report the expressions of the d(q) term: 

d{q)=dTr[(i + l> )l x hhls]- 

No other g-dependent structure can appear at the level of 4-point functions. For the 3-point 
functions, there are six possible terms that can contribute to c(q): 

Jmax 

3=1 



In the renormalizable gauge, j max = 3. The 2- and 1-point spurious functions b(q) and d(q) 
contain eight and four terms, respectively 

After fixing the form of all the spurious terms, our calculation is reduced to the algebraic 
problem of extracting all the coefficients. This is achieved simply by evaluating N(q) at different 
values of the integration momentum q. As a further simplification, there is a very good set of 
such points: we can use values of q for which a subset of denominators Di vanish. Operating 
in this manner, the system becomes "triangular" ; we can solve first for 4-point functions, then 
3-point functions and so on. 

To conclude this section, let us summarize the recipe for the calculation. We first calculate 
all the coefficients in Eq. ([2]), by evaluating the numerator of the integrand N(q) for a set 
of values of the integration momentum q. Note that we do not need to repeat this for all 
Feynman diagrams: we can group them and expand for (sub)amplitudes directly. We just need 
to specify external momenta, polarization vectors and masses and proceed with the reduction. 
Concerning N(q) we can choose how to proceed according to the specific calculation: as an 
interesting development, N(q) could be determined numerically via recursion relations. Once 
all the coefficients have been determined, we can multiply them by the corresponding scalar 
integrals. For massive integrals we use FF by G. J. van Oldenborgh [14J, while in the massless 
case the integrals are evaluated using the code OneLOop by A. van Hameren |15j . 



3 Computation of rational terms 

As mentioned in the previous section, the reduction has been performed assuming a purely four 
dimensional numerator (this singles out the so called cut-constructable part of the amplitude) . 
In this section we describe our method to calculate the rational parts of the amplitude. In our 
approach, these contributions originate from the fact that up to now we expressed the numerator 
in terms of Di, while the functions appearing in the denominator are n-dimensional Di. Let us 
go back to the integrand A{q) of Eq. (JTJ) and insert the expression for N(q) of Eq. ([2]), that we 
obtained after determining all the coefficients for both for regular and spurious terms. Now, 
rewriting Di by means of 
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we obtain 
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The rational part is produced, after integrating over cPg, by the q 2 dependence in Z{. The 
expressions for all relevant integrals are reported in the Appendix of Ref. [12J. 

In addition to the one just described, there might be other sources of rational terms coming 
from objects of dimension e in the numerator N(q), according to the specific calculation. These, 
for example, could originate from the contraction of Dirac matrices or from powers of q 2 [16]. 



4 Numerical Tests 



As an example of applications of the method, we calculated 4-photon and 6-photon amplitudes, 
via fermionic loop of mass m/ |12j . Those processes do not have an immediate physical interest, 
however they provide a playground for testing and comparing different methods of calculation. 

In order to perform the calculation, for each phase-space point we should provide as input 
parameters: the external momenta pi (for this particular example they are massless, i.e. pf = 0); 
the masses of propagators in the loop (fermion masses mj in this example); the polarization 
vectors for the photons: here we use various helicity configurations. 

For the 4-photon amplitudes, we performed a comparison with the analytic result presented 

by Gounaris et al. [T7J. We compared separately the rational part and the cut-constructable 

part, both in the massive and massless cases, finding perfect agreement. We repeated the 

f f f 

calculation for the three helicity configurations F^_ +++ , Fi,, and F{, presented in the 

Appendix of Ref. [IT] . 

Coming to the 6-photons, there are few results available in the literature for the massless 
case. Some time ago, Mahlon presented an exact analytic result for the helicity configuration 

[+ H ] [18] . More recently, Nagy and Soper, with a fully numerical approach, obtained 

the results for the configurations [+ H ] and [+ H 1 — ] [19]. The same results were 

also recently presented by Binoth et al. [20], that also provide analytic expressions. Our results 
are in full agreement with all previous calculations. We also checked that the cut-constructable 

part for configurations [H ] and [ ] and the rational terms for all helicity 

configurations are identically zero. 

Finally we calculated the 6-photon amplitudes with massive internal fermions. Also in this 
case there is no contribution from rational parts. The results have been presented in Ref. |12j . 



5 Summary and Conclusions 

The discovery potential of LHC requires NLO calculations. At present, there is a variety of 
interesting options available to perform one-loop multi-leg calculations, however not a universal 
method. 

We recently proposed a method for the numerical evaluation of one-loop amplitudes in 
which the reduction is performed at the integrand level. The method is rather "young" and 



there is still room for improvements and optimizations, that are currently in progress. The first 
results are encouraging. We plan to present more results soon, for processes of interest at future 
experiments. 
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